Speckle visibility spectroscopy and variable granular fluidization 
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We introduce a dynamic light scattering technique capable of resolving motion that changes 
systematically, and rapidly, with time. It is based on the visibility of a speckle pattern for a given 
exposure duration. Applying this to a vibrated layer of glass beads, we measure the granular 
temperature and its variation with phase in the oscillation cycle. We observe several transitions 
involving jammed states, where the grains are at rest during some portion of the cycle. We also 
observe a two-step decay of the temperature on approach to jamming. 
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A wealth of spectacular phenomena occur when gran- 
ular materials are subjected to periodic vertical vibra- 
tion jl], ||. For example, shallow layers exhibit pe- 
riod doubling and pattern formation; deep layers ex- 
hibit heaping, convection, and logarithmically-slow com- 
paction. It is relevant to ask: What grain-scale physics 
gives rise to all this intriguing macroscopic behavior? 
Since the driving is generally at high-amplitude but low- 
frequency, the dynamics vary dramatically during the os- 
cillation cycle. At first, the grains are jammed || 
completely at rest in some random packing configura- 
tion. When the downward acceleration exceeds —g, the 
grains are launched upward from the plate. The layer 
expands as the grains collide and move about randomly. 
Soon they crash back into the bottom plate, and pick 
up energy from this impact. Finally, they come to rest 
after rattling away energy by inelastic collisions. This 
approach to jamming can be accompanied by clustering 
and "inelastic collapse" , a finite time singularity where 
the collision rate diverges and the collision length van- 
ishes H |(| . This regime is especially significant because 
granular hydrodynamics |?], [|, ^ and statistical mechan- 
ics, extended to an athermal system, both break down. 

Unfortunately, the above sequence of dynamics has not 
been experimentally accessible. Because of multiple light 
scattering, video imaging is restricted to dilute granu- 
lar gasses or surface behavior (To). Even then, the spa- 
tial resolution is much larger than the collision length as 
the grains come to rest. This limitation holds for other 
imaging techniques as well, like MRI ||], x-ray mi- 
crotomography [ fl3| , and positron emission particle track- 
ing |Q. Furthermore, none is fast enough to capture the 
high collision rates when the grains are barely fluidized. 
By contrast, diffusing-wave spectroscopy (DWS) p5) is 
a dynamic light scattering (DLS) jl6| method that ap- 
plies to bulk granular media. It has superior spatial and 
temporal resolution, and can be extended to unsteady 
dynamics [ p7| . However, it is based on temporal correla- 
tion functions, which implicitly assume that all times are 
statistically equivalent; therefore, it is not appropriate for 
periodic or aging systems. Altogether, the leading probe 



of dynamics in vibrated 3D granular systems currently 
is an NMR technique [jl2|. It allows tracking of individ- 
ual grains throughout a highly-fluidized sample to within 
about 150 and 1.4 ms. 

In this paper, we introduce a new DLS technique and 
use it to study the bulk grain dynamics throughout the 
oscillation cycle. Here the resolution is limited by the 
wavelength of light and by the speed of a fast CCD cam- 
era. Taking advantage of multiple light scattering, we 
achieve a resolution of «1 nm and «20 ^s. With this ad- 
vance, we have unlocked regimes where the grains barely 
move. In particular, we observe three dynamic transi- 
tions: the onset of fluidization, where the acceleration 
amplitude just exceeds g; a jamming transition, where 
the grains crash into the plate and then come to rest; and 
a transition to continuous fluidization, where the grains 
do not jam up at any point during the cycle. In contrast, 
the NMR study of Ref . Jl2| was conducted far above this 
point, where granular hydrodynamics is applicable in the 
bulk. Our three transitions all involve a jammed state, 
and therefore cannot be captured by granular hydrody- 
namics. Since our observations quantify the microscop- 
ies that underlie a host of intriguing phenomena, they 
present a theoretical challenge. 

Our granular system consists of 780 ± 35 /im diam- 
eter glass beads, approximately twelve layers deep in a 
10 x 10 cm 2 box with a flat transparent bottom, ver- 
tical walls, and open top. This is mounted on a shake 
table, which in turn is leveled on an optical bench. Two 
three-axis accelerometers monitor the quality and peak 
amplitude, a = Tg, of the vertical sinusoidal oscilla- 
tions. All data are taken at frequency / = 10 Hz. We 
define the phase to be <f> — when the plate is at height 
z = and moving upward. Properly leveled, we observe 
no discernable heaping, pattern formation, or convection 
across the range of amplitudes studied, < V < 2.2. 

To measure grain motion, we introduce a method that 
we call "speckle visibility spectroscopy" (SVS). As shown 
in Fig. 0a, we illuminate the grains from above with a 
wl cm diameter beam of a 100 mW frequency-doubled 
Nd-YAG laser (A = 532 nm) . Photons perform a random 
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walk with transport mean free path of «4 grains. Since 
the sample is 12 layers deep, about ^ of the photons reach 
the bottom plate, and the rest are backscattered after a 
few scattering events. The backscattered light forms a 
speckle pattern in the far field, which we detect using a 
digital linescan CCD camera (Basler-160: 1024 pixels, 8 
bits deep), and no intervening optics except a 532 nm 
filter. The sample-CCD distance is «20 cm, such that 
the speckle size is comparable to the pixel size (10 fixn). 
Therefore, as the grains move and the speckle pattern 
changes, large intensity fluctuations occur at each pixel. 
For the speckle pattern to be visible, however, the ex- 
posure time of the CCD must be short compared to the 
time scale for speckle fluctuations. If the exposure time 
is longer, then the speckle blurs out and the same average 
intensity is recorded at each pixel. This is the essence of 
SVS. It is illustrated qualitatively in Fig. 0b by the space- 
time plot of speckle vs phase during the cycle. When the 
grains are at rest on the plate, the speckle is clearly visi- 
ble. When the grains are fluidized, the speckle is blurry. 
The speckle is least visible just after impact, where rapid 
grain motion is excited by the sudden injection of en- 
ergy. This is reminiscent of "laser speckle photography" , 
where the absence of speckle in a laser-illuminated scene 
indicates motion |l8|, |l9[ ]. 

The key measurable quantity in SVS is the variance of 
intensity across the pixels. For an exposure duration T, 
each of the N pixels reports a time-integrated intensity, 
Si = Jo h(t)dt/T. The ensemble -averaged intensity and 
intensity-squared are computed as (I)t — J^iLi^i/N 

and (I 2 )t = J2iLi Si 2 /N, respectively. If there are 
enough speckles, then the former is independent of T 
and the subscript may be dropped. By contrast, (I 2 )t 
depends on T and indicates the visibility of the speckles 
to the extent that it exceeds (I) 2 - To quantify visibility 
on a scale of — 1, we define a normalized variance: 



V 2 (T)^[(I 2 ) T /(I) 2 -l] 



(1) 



The factor 1/(3 is roughly the number of speckles per 
pixel, and is determined experimentally by measuring the 
system at rest. Fig. |l|c shows an example of Vz(T) vs 
phase during the cycle, for a fixed exposure of T = 100 fis. 
It is closer to when the grains are moving rapidly, and 
closer to 1 when the grains are coming to rest. 

To relate the variance to grain motion, note that 
Eq. ([!]) involves ensemble averages, rather than time av- 
erages. Therefore, the Siegert relation |ll| holds, giving 
(I 2 ) T = <f T / T I % {t')h(t")dt'dt"/T 2 ) t = (J> a /o r J> + 
P\gi(t' - t")\ 2 dt'dt" /T 2 . Here, f3 is the same as above, 
and gi(t) is the normalized electric field autocorrelation. 
Since g\ (t) is even, the double integral simplifies and we 
arrive at the fundamental equation of SVS: 
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V 2 (T) = 2 [ T (l-t/T)\ 9l (t)\ 2 dt/T. 
Jo 
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FIG. 1: (a) Theoretical position of the cell, bead layer, and 
their difference (5x), vs phase <j> for a sinusoidally oscillated 
cell. The particular curves are for V — 1.25 and / = 10 Hz. 
For intuition, the bead layer is modeled as a slab with an 
effective coefficient of restitution; a value of e = 0.5 is used for 
the data shown. The inset shows the optical geometry, (b) A 
representative speckle image vs <j> gathered by the linear CCD 
for T = 1.25, / = 10 Hz, and T = 100 /is. For clarity, only 
200 of 1024 pixels are shown, (c) The averaged normalized 
variance V2 of 300 full images versus <j> for the above settings. 



The variance is thus a weighted average of |ffi(i)| 2 over 
the exposure time T. Both functions are 1 at short times, 
and at long times. Given gi(t) from SVS measure- 
ments, the scattering site motion may then be deduced 
by standard DLS practice. For random ballistic motion 
of average speed 6v, the theory of DWS for backscattered 
light gives g\{t) — exp(— 7*) with 7 = inSv/X |15| ]. The 
corresponding variance is 

V 2 (T) = 2[exp(-2 7 T) - (1 - 2 7 T)]/(27T) 2 . (3) 

At short times, the initial decay is linear: V 2 (T) 1 — 
37T; at long times, it is a power-law: V 2 (T) w 1/(jT). 
The heavy weighting in Eq. (Q) near t = slows the 
decay, aiding in the measurement of fast processes. 
We now return to Fig. |l|, where the grains are vibrated 
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FIG. 2: (a) Semilog plot of Vi vs exposure time T for F = 1.25 
and / = 10 Hz at three phases (j>. The symbols correspond to 
the markers in Fig. lc. The curves are single parameter fits 
using Eq. (||); only the data denoted by the larger symbols 
has been fit. The data fit by the dashed curve represent the 
effect of speckle wash alone, (b) The fitted rates vs phase. 



at / = 10 Hz and T — 1.25, and exploit our SVS method. 
Three representative phases in the cycle are marked by 
triangles in the 1^(100 fis) data in Fig. ^jc: where the 
grain motion is most rapid both in mid-flight and after 
impact, and where the grain motion ceases. For each of 
these three events, we show ^(T) vs T in Fig. ||a. All 
data tend toward 1 (0) at short (long) times, as expected, 
where the speckle is most (least) visible. The actual dy- 
namic range of our data is limited by two effects that are 
specific to our experiment (not the SVS method). First, 
for exposures faster than w 50 /is, our 100 mW laser 
produces a signal that is a small fraction of the 0-255 
range of the CCD. Therefore, the distribution of intensi- 
ties is binned coarsely, which systematically distorts the 
variance. This limitation could be reduced by a brighter 
laser. Second, the macroscopic motion of the system con- 
tributes to the variance, becoming significant for T —> oo 
and/or Vi — > 1. In particular, even if the grains are at 
rest, the speckles form a static pattern that washes as a 
whole across the CCD pixels. This "speckle-wash" can 
be seen as a swirling in the space-time plot of Fig. |l|b. 

The reliable portion of the variance data in Fig. |2|a 
compares well with Eq. (0) , as shown by the curves. The 
form of the decay is therefore consistent with expectation: 
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FIG. 3: Grayscale plot of decay rate 7, based on fits to Eq. (Q) , 
as a function of both peak acceleration F and phase (j> in the 
oscillation cycle; black=0, white=120 kHz. The solid white 
curve indicates where the grains are launched, at an instan- 
taneous acceleration of — g. The dotted white curve indicates 
when a free-fall object would land. 



the grains appear to undergo random ballistic motion 
with some average fluctuation speed Sv. For granular 
materials, the kinetic energy associated with Sv is called 
the granular temperature; it is directly proportional to 
the decay rate 7 = AttSv/X. Results for 7 are shown in 
Fig. |^b as a function of phase in the cycle. The grains 
"heat up" both in mid- flight and from impact. In both 
cases, they lose energy through inelastic collisions. About 
■i cycle (20 ms) after impact, the grains lose their energy 
and come to rest in a jammed state. Due to speckle 
wash, however, 7 doesn't decay to zero but rather to a 
readily-identifiable baseline. To our knowledge, this is 
the first measurement of the bulk granular temperature 
during the collapse. 

We now repeat the experiment versus T, at constant 
/ = 10 Hz. A grayscale plot of 7 is shown in Fig. || 
vs both r and <f). Effectively, this is a phase diagram 
denoting the relative fluidization of the medium at dif- 
ferent forcing rates and at different points in the cycle. 
Several transitions can be observed. First, below T = 1, 
there is no fluidization at all. Above r = 1, the grains 
become fluidized for some portion of the cycle. In par- 
ticular, the grains are launched from the plate when the 
instantaneous acceleration is —g, as denoted by a solid- 
white curve. The dotted-white curve denotes where the 
grains crash back down, assuming they undergo simple 
free-fall. Some time after the grains land, they lose their 
energy via collisions and jam. The onset of both jam- 
ming and unjamming depend on the driving amplitude. 
For r > 1.6, they merge together and the system un- 
dergoes a transition to continuous fluidization. Though 
the grains never come to rest, their dynamics still vary 
throughout the cycle. As expected, the motion is fastest 
(slowest) just after (before) impact. 
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FIG. 4: Average fluctuation speed vs time after impact, ex- 
tracted from Fig. ^, for different accelerations T. For the solid 
curves, all below the transition to continuous fluidization, 5v 
increases monotonically with driving amplitude. Above the 
transition, Sv never goes to zero and doesn't spike as strongly 
when the grains impact the plate. 

Both the process of jamming, and the nature of the 
transition to continuous fluidization, may be studied in 
terms of the fluctuation speed, 8v, vs the time following 
impact, At. These data are extracted from Fig. || and dis- 
played in Fig. ||. For amplitudes 1 < T < 1.6, 6v spikes 
at impact and then decays to a plateau after «2 ms; the 
initial decay rate appears to be independent of T. Fol- 
lowing the plateau, 5v decays to zero (within resolution 
set by speckle wash) as the grains jam up. This second 
decay is slower than the first. The temperature spike, 
the level of the plateau, and the duration of the plateau, 
all increase monotonically with T. The impacts become 
progressively more violent, until suddenly at T > 1.6 the 
plateau extends across the entire cycle and the grains 
never come to rest. The impacts are no longer as vio- 
lent; instead, the grains stay fluidized and Sv exhibits a 
smoother variation with time. 

Conclusion The use of area detectors for "multi- 
speckle" dynamic light scattering is on the rise [^0[ |2l], 
|22| , |2l| . The general approach has been to autocorrelate 
each pixel, and then to average the results afterwards, 
all in software. Our new method, speckle variance spec- 
troscopy (SVS), is dramatically different and offers ad- 
vantages in terms of both simplicity and applicability. In 
effect we have created a "speckle ensemble correlator" 
in which all averages are explicitly-computed ensemble 
averages of a single exposure of the CCD camera; nei- 
ther time averages, temporal autocorrelations, nor image 
storage are necessary. While the time resolution of prior 
methods is currently no better than s=s 2 ms, we achieve 
~ 20 /is. More significantly, we can follow dynamics 
that change on equally rapid time scales. This unprece- 



dented resolution allows us to capture the transition to 
continuous fluidization in Fig. ^ and the two-step decay 
of the granular temperature in Fig. ||. With a brighter 
laser, and better control of speckle wash, we are poised 
to study these and other granular dynamics in greater 
detail. SVS also opens a new window for the study of 
bubble rearrangements in coarsening foams, motion and 
aging in glassy suspensions, gelation, phase-separation, 
and other phenomena that exhibit fast or quickly evolv- 
ing nonstationary dynamics. 
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